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Abstract 

We  demonstrate  that  high-frequency  cooperative  dynamics  of  an  array  of  coupled  nonlinear  spin- 
torque  nano-oscillators  (STNO)  can  be  controlled  by  introduction  of  an  additional  external  phase 
shift  /?c  between  the  microwave  current,  which  couples  STNOs,  and  the  total  array  microwave 
voltage.  When  this  external  phase  shift  /3C  compensates  the  intrinsic  phase  shift  /?o,  caused  by  the 
STNO  nonlinearity,  a  phase-locking  regime  with  increased  output  power  and  vanishing  inhomoge¬ 
neous  linewidth  broadening  is  achieved.  In  the  opposite  case,  when  external  and  intrinsic  phase 
shifts  are  added,  the  STNO  array  demonstrates  a  frustration  regime  with  low  output  power  and 
wide  and  noisy  frequency  spectrum. 

PACS  numbers:  85.75.-d,  75.30.Ds,  75.40.Gb 
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The  major  problems  in  the  practical  realization  of  spin-torque  nano-oscillators  (STNO) 
[1-4]  are  their  low  output  power  and  large  generation  lincwidth.  Both  these  problems  can 
be  solved  by  using  the  phenomenon  of  mutual  phase- locking  in  an  array  of  coupled  STNOs 
[5-9], 

In  contrast  with  the  majority  of  conventional  oscillators,  STNOs  are  strongly  nonlinear: 
both  the  generated  frequency  u>(p)  and  the  total  damping  T(p)  strongly  depend  on  the  gen¬ 
erated  power  p  (see  [10]  for  details).  It  is  convenient  to  characterize  the  STNO  nonlinearity 
by  a  dimensionless  parameter  (see  Eq.  (33)  in  [10]): 

duj(p)/dp 
dT  ( p )  j  dp  ’ 

which  for  typical  STNO  parameters  could  have  the  magnitude  as  large  as  \v\  ~  100. 

The  nonlinearity  significantly  influences  the  interaction  of  STNOs  with  external  signals 
and  with  each  other.  In  particular,  it  substantially  increases  the  frequency  band  A  of 
phase-locking  of  STNO  to  an  external  signal  (compared  to  the  linear  value  Ao  of  this  band) 
[8,  11,  12]  and  it  creates  a  large  (about  ±90°)  intrinsic  phase  shift  /30  [13-15]  between  the 
driving  signal  and  the  oscillations  of  STNO  resistance  (see  Eqs.  (52),  (53)  in  [10]): 

A  =  Vl  +  zz2A0  ,  (30  =  —  arctan(i/)  .  (2) 

The  intrinsic  phase  shift  (do  does  not  change  the  conditions  of  phase-locking  of  an  STNO 
to  an  external  signal,  the  phase  of  which  is  independent  of  the  phase  of  STNO,  but  may 
be  of  a  critical  importance  for  the  process  of  mutual  phase-locking  in  coupled  STNO  array, 
where  the  driving  signal  is  created  by  all  the  STNOs  in  the  array. 

In  this  work  we  demonstrate  that  the  large  intrinsic  phase  shift  (do  can  completely  suppress 
mutual  phase- locking  in  an  STNO  array.  We  also  show  that  the  disruptive  influence  of  (do 
can  be  compensated  by  introduction  of  an  additional  external  phase  shift  (dc  of  a  proper  sign 
between  the  total  microwave  current  driving  the  array  (and  created  by  all  the  STNOs  in  the 
array)  and  the  total  microwave  voltage.  The  introduction  of  this  additional  phase  shift  (dc 
does  not  affect  the  nonlinear  enhancement  of  the  phase-locking  band  (see  first  equation  in 
(2))- 

The  mutual  phase-locking  of  STNOs  was  studied  numerically  on  the  example  of  an  array 
shown  in  Fig.  1.  It  consists  of  N  =  10  nano-pillar  STNOs  electrically  connected  in  series, 
biased  by  the  dc  current  J0  =  1  mA,  and  connected  in  parallel  with  an  external  RLC  circuit. 
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The  dynamics  of  each  STNO  is  described  by  the  Landau-Lifshits-Gilbert-Slonczewski 
equation  [1]  for  the  magnetization  vector  M,  of  ’’free”  layer  of  i-th  STNO: 

Mi  =  7[. Hi  x  Mi]  +  j^r[Mi  x  Mi]  (3) 

Here  7  =  2n  ■  2.8  MHz/Oe,  HL  =  (Hq  —  4n MitZ)z,  Hq  =  10  kOe  is  the  bias  magnetic  field 
directed  perpendicularly  to  the  ’’free”  layer  (along  coordinate  z),  and  term  — 47T MitZ  is  the 
dipolar  field  responsible  for  the  nonlinearity  of  STNO.  The  second  term  in  Eq.  (3)  is  the 
Gilbert  damping  term,  «g  =  0.01  and  47 rM0  =  8  kOe  are  the  Gilbert  damping  constant  and 
saturation  magnetization  of  the  ’’free”  layer,  respectively. 

The  last  term  in  Eq.  (3)  is  the  spin-transfer  torque  which  describes  the  action  of  the 
electric  current  /(f)  on  the  i-th  STNO.  Here  <7*  =  £i'yh/(2eM0Vo),  £t  is  the  dimensionless 
spin-polarization  efficiency  for  i-th  STNO,  h  =  h/( 2n)  is  the  reduced  Planck  constant,  e  is 
the  modulus  of  electron  charge,  Vo  =  2.35  x  104  nm3  is  the  volume  of  the  ’’free”  STNO  layer 
(which  corresponds  to  the  circular  nano-pillar  of  the  radius  50  nm  and  thickness  3  nm),  and 
p  =  sin  dp  x  +  cos  6pz  is  the  direction  of  spin  polarization  of  the  current  (9P  =  30°  in  our 
calculations). 

Spin-polarization  efficiencies  £i  were  chosen  to  be  different  for  different  STNOs  (1  % 
spread  around  the  central  value  of  £0  =  0-4),  thus  creating  a  spread  of  the  free-running 
frequenies  of  uncoupled  STNOs  in  the  interval  A/0  ~  100  MHz  around  the  central  generation 
frequency  of  /0  ~  12.16  GHz  (see  Fig.  2a  and  discussion  below). 

The  electric  current  /(f)  consists  of  the  dc  bias  current  Iq  and  microwave  current  A /(f) 
(see  Fig.  1)  created  by  the  magnetization  precession  of  each  STNO  due  to  giant  magnetore¬ 
sistance  (GMR)  or  tunneling  magnetoresistance  (TMR)  effect.  For  the  considered  electrical 
scheme  Fig.  1  A /  =  dQ/d,t,  where  Q  is  the  charge  on  the  capacitance  C,  which  obeys  the 
equation 

LQ+  R  +  Q  +  Q/C=-J2Ri( f)/o,  (4) 

i  J  i 

where  Ri(t)  is  the  resistance  of  the  i-th  STNO: 

Ri(t)  —  Ro  +  A  Riit)  =  Rq  1  +  (1  —  cos  #j(t))  .  (5) 

Here  R$  =  10  Ohm  is  the  electric  resistance  of  STNO  in  parallel  magnetic  state,  /x  =  0.3 
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is  the  magnetoresistance  of  the  oscillators,  and  6{(t)  =  arccos(Af?:(t)  ■  p/Mq)  is  the  time- 
dependent  angle  between  the  magnetization  directions  in  the  ’’free”  and  ’’fixed”  magnetic 
layers. 

The  microwave  current  A /(f),  created  by  all  the  STNOs  and  acting  on  each  of  them, 
serves  not  only  as  the  output  signal  from  the  STNO  array,  but,  also,  is  the  coupling  signal 
that  can  drive  the  array  of  oscillators. 

The  role  of  the  RLC  circuit,  placed  in  parallel  with  the  STNO  array,  was  to  create 
an  external  phase  shift  (3C  =  arcta n((/c  —  /0)/A/c)  (where  fc  =  1/(27 Ty/LC),  A fc  = 
(R  +  NR0)/(4:'jtL)  =  1  GHz,  and  f0  is  the  generation  frequency  of  the  array)  between  the 
microwave  current  A I(t)  flowing  through  the  array  and  the  generated  microwave  voltage 

e{t)  =  ioEi*Ri(t)- 

The  results  of  the  numerical  modeling  of  the  dynamics  of  the  STNO  array  Fig.  1  are 
presented  in  Fig.  2,  which  demonstrates  temporal  profiles  of  the  microwave  current  A I(t) 
(left  column)  and  spectral  density  of  power  on  the  active  load  R  =  50  Ohm  (right  column) 
for  several  different  parameters  of  the  RLC  circuit. 

The  upper  row  (Fig.  2(a))  shows  the  case  of  uncoupled  STNOs,  when  the  feedback  action 
of  the  common  microwave  current  A /(f)  on  the  STNO  dynamics  was  neglected.  The  power 
spectrum  in  this  case  consists  of  N  =  10  narrow  peaks  generated  by  each  of  STNOs.  The 
incoherent  sum  of  N  =  10  independent  oscillations  with  different  frequencies  leads  to  chaotic- 
like  temporal  profile  of  the  current  A I(t)  (see  Fig.  2(a),  left  column). 

Fig.  2(b)  shows  the  case  of  coupled  STNO  array  when  the  resonance  frequency  of  the 
RLC  circuit  was  exactly  equal  to  the  central  generation  frequency  fc  =  fo  =  12.16  GHz 
(the  external  phase  shift  was  fic  —  0). 

It  is  important  to  stress,  that,  if  one  ignores  the  intrinsic  phase  shift  flo,  the  generated 
current  A I  ~  20  /i  A  is  sufficient  to  have  phase-locking  in  the  frequency  interval  ~  300  MHz 
[8],  which  is  several  times  larger  than  the  spreading  of  the  free-running  frequencies  of  our 
STNO  array.  Thus,  naively,  one  may  expect  that  the  coupled  STNO  array  will  be  completely 
phase-locked. 

The  numerical  results  Fig.  2(b),  that  implicitly  take  into  account  the  phase  shift  /So,  prove 
otherwise.  The  temporal  profile  of  the  current  A /(f)  is  still  quasi-chaotic  (left  column),  while 
its  power  spectrum  is  wide  and  noisy,  and  is  shifted  to  the  region  of  lower  frequencues  (right 
column).  Similar  shift  (only  to  the  region  of  higher  frequencies)  was  found  in  simulations 
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[16]  performed  for  the  case  of  in-plane  magnetized  STNO  array  (i.e.,  STNOs  having  opposite 
sign  of  the  nonlinearity  v  Eq.  (1)). 

The  reason  for  this  incoherent  behavior  of  the  array  is  the  large  intrinsic  phase  shift  /do, 
that  for  the  parameters  of  our  array  can  be  evaluated  using  [10]  as  /3q  ~  —  arctan(z/)  « 
—89.4°.  Due  to  this  large  phase  shift  the  oscillations  in  each  STNO  dynamically  adjust  their 
phases  roughly  90°  behind  the  phase  of  the  current  A /(f),  and,  since  the  driving  current 
A I  ft)  itself  is  created  by  the  STNOs,  such  a  90°  phase  delay  destroys  any  coherence  that 
can  spontaneously  arise  in  the  array. 

A  very  different  picture  is  seen  when  one  introduces  an  external  positive  phase  shift 
/3C  =  arctan((/c  —  /0)/A/c)  =  45°,  that  partly  compensates  the  intrinsic  phase  shift  /30  to 
give  the  total  phase  shift  of  (3  =  (3o  +  /3C  ~  —45°.  In  that  case  (see  Fig.  2(c))  the  array 
enters  a  coherent  phase-locked  generation  regime  with  narrow  frequency  spectrum  and  large 
output  power. 

To  demonstrate  that  the  observed  qualitative  change  of  the  STNO  array  dynamics  is 
caused  solely  by  the  phase  mechanism,  we  performed  simulations  for  the  case  when  the 
external  phase  shift  introduced  by  the  RLC  circuit  is  negative  /3C  =  —45°,  thus  making  the 
total  phase  shift  /3  =  /3o  +  /3c  ~  —135°  (see  Fig.  2(d)).  In  this  case  the  dynamics  of  the  STNO 
array  qualitatively  differs  from  all  the  previously  considered  cases.  The  generation  spectrum 
remains  chaotic  (as  in  Fig.  2(b)),  but  the  amplitude  of  the  current  A /  is  drastically  reduced. 
This  reduction  is  not  related  to  the  increase  of  the  RLC  circuit  impedance  (which  was  exactly 
the  same  in  the  cases  (c)  and  (d)).  Our  additional  simulations  with  a  different  number  N 
of  oscillators  indicate  that  the  magnitude  of  the  current  decreases  with  the  increase  of  the 
number  of  interacting  oscillators  and  completely  vanishes  in  the  limit  N  — >  oo. 

This  counter-intuitive  result  can  be  understood  by  noting  that  for  the  large  phase  shifts 
|/3 1  >  90°  the  interaction  between  the  STNOs  has  a  repulsive  nature,  i.e.  oscillators  dynam¬ 
ically  adjust  their  phases  to  minimize  the  mutual  coupling.  Therefore,  the  observed  drastic 
reduction  of  the  coupling  current  A /  in  this  frustration  cooperative  regime  is  a  result  of  the 
dynamical  destructive  interference  of  oscillations  of  different  STNOs. 

Our  simulation  results  can  be  qualitatively  explained  using  the  modified  Kuramoto  model 
of  coupled  auto-oscillators,  which  takes  into  account  phase  shift  /3.  In  the  framework  of  this 
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model,  each  STNO  is  characterized  by  only  its  phase  (f>i(t)  (see  [10]  for  details): 


4>i  =  uji  +  A  ^2  sin(0j  ~<f>i  +  P)  ■  (6) 

3 

Here  u>i  =  27r/,;  are  the  free-running  oscillator  frequencies,  A  >  0  describes  the  interaction 
strength  between  the  oscillators,  and  /3  =  /30  +  /3C  is  the  total  phase  shift. 

In  the  limit  of  large  coupling  A  — >■  oo  the  system  (6)  has  the  coherent  stationary  solution 

4>iit)  —  ut ,  uj  ~  luo  +  AT  A  sin  /3  .  (7) 

Eq.  (7)  explains  the  above  mentioned  shift  of  averaged  frequency  u  in  the  coupled  STNO 
array.  Linear  stability  analysis  shows  that  the  coherent  solution  Eq.  (7)  is  stable  only  for 
cos/3  >  0,  i.e.  for  \f3\  <  90°.  For  larger  angles  |/3|  >  90°  (cos/3  <  0)  the  coherent  dynamics 
of  STNO  array  is  impossible  independently  of  the  interaction  strength  A.  In  this  case  the 
incoherent  (or  frustrated )  state  with  evenly  distributed  phases  </>*  becomes  stable.  The  total 
output  power  (oc  |  exp[i(/>j(f)]|1 2 3 4 5),  vanishes  completely  in  this  frustrated  state. 

In  conclusion,  we  have  shown  that  the  cooperative  dynamics  in  an  array  of  coupled 
STNOs  can  be  controlled  by  the  introduction  of  an  additional  external  phase  shift  /3C  to  the 
microwave  coupling  signal.  By  changing  /3C,  STNO  array  can  be  driven  into  one  of  three 
qualitatively  different  regimes:  coherent  phase-locked  regime  (see  Fig.  2(c))  with  large  output 
power  and  narrow  generation  lincwidth,  intermediate  chaotic  regime  (see  Fig.  2(b))  with 
wide  spectrum  and  relatively  large  output  power,  and  the  frustration  regime  (see  Fig.  2(c)) 
with  a  vanishingly  small  output  power.  We  believe  that  our  results  will  be  important  for 
the  design  of  practical  oscillators  based  on  the  phase-locked  arrays  of  STNOs. 
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L564),  and  from  the  U.S.  Army  Research  Office  (MURI  grant  No.  W911NF-04-1-0247). 
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Figure  captions 


FIG.  1:  Scheme  of  the  studied  STNO  array  consisting  of  N  oscillators  connected  in  series  and 
biased  by  the  dc  current  Iq.  The  microwave  current  A I(t),  created  by  the  array,  flows  through  the 
external  LCR  circuit.  By  varying  the  capacitance  C  of  this  circuit  one  can  control  the  phase  shift 
j3c  between  the  microwave  current  A I(t)  and  the  total  microwave  voltage  £(t)  =  Io  Yli  ARi{t)  on 
the  array. 


FIG.  2:  Typical  temporal  profiles  of  the  microwave  current  A I(t)  (left  column)  and  power  spectral 
density  on  the  active  load  R  (right  column)  for:  (a)  -  uncoupled  STNOs  and  (b),  (c),  (d)  - 
coupled  STNO  array  with  different  resonance  frequencies  fc  of  the  RLC  circuit:  (b)  fc  =  fo,  which 
corresponds  to  the  zero  external  phase  shift  (3C  =  0;  (c)  fc  —  fo  +  Afc,  ( positive  external  phase 
shift  (3C  =  45°);  (d)  fc  =  fo  ~  Afc,  ( negative  external  phase  shift  f3c  =  —45°).  Blue  lines  in  the  left 
panel  show  the  envelopes  of  the  microwave  current. 
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